Perturbative treatment of inter-site couplings in the local descrip¬ 
tion of open quantum networks 


A. S. Trushechkin and I. V. Volovich 

Steklov Mathematical Institute of the Russian Academy of Sciences - Gubkina 8, Moscow 119991, Russia 







PACS 03.65. Yz - Decoherence; open systems; quantum statistical methods 
PACS 05.60. Gg - Quantum transport 

Abstract - The problem of construction of a quantum master equation for a system of sites 
weakly coupled to each other and to one or more reservoirs (open quantum network) is considered. 
Microscopic derivation of a quantum master equation requires a diagonalization of the Hamiltonian 
of the network, which can be a difficult task. When the inter-site couplings are weak, the local 
approach, which neglects the influence of the inter-site couplings on the system-reservoir couplings, 
is often used. Recently, some doubts were cast to the consistency of the local approach. We develop 
a systematic perturbation expansion to derive corrections to the local approach and establish its 
range of validity. Using this extension of the local approach, we derive an expression for the heat 
flux for a particular model and show that it does not violate the second law of thermodynamics. 
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Introduction. — Quantum master equation frame¬ 
work is commonly used in the description of charge and 
energy transfer in quantum networks [1] . Quantum master 
equation can be microscopically derived, but this requires 
a diagonalization of the Hamiltonian of the network. This 
can be a difficult task. Even if not, this master equation 
can be complicated and, hence, difficult for analysis. This 
approach is sometimes called global, because it considers 
the network as a whole. When the inter-site couplings 
in the network are weak, it seems reasonable to neglect 
their influence on the system-reservoir couplings. This 
approach is sometimes called local and is often used as a 
simplification, since it does not require a diagonalization 
of the full network Hamiltonian. 


However, recently several authors US] reported the dis¬ 
crepancies between the local and global descriptions of 
open quantum systems and, in particular, the violation of 
the second law of thermodynamics by the local approach 

m- 


We show that the viewpoint of perturbation theory al¬ 
lows to justify this approach, establish its range of validity 
and derive corrections to it. The small parameter is the 
inter-site coupling constant. In this formalism, we derive 
an expression for the heat flux for the model analysed in 
[2 and show that it does not violate the second law of 
thermo dynamics. 


Global and local quantum master equations. — 

Consider a large system consisting of a relevant system 
(“system”) and a reservoir. The Hamiltonian has the form 

i? = iIs + iJjj + Ai?7. (1) 

Here Ms and Mr specify the free dynamics of the system 
and the reservoir accordingly. The interaction is speci¬ 
fied by the term Mj = where Aa act on 

the system, Ma act on the reservoir, and A is a small di¬ 
mensionless parameter. The application of weak coupling 
limit with respect to A leads to the following master equa¬ 
tion for the reduced density matrix p = p[t) of the system 

[5H2]: 

^^=-^[Ms + X^MLs.p\ + >?V{p), 

where 

'^(p) = T.T. 'Jaa' (^) (^)'^ 

uj a,ex' 

-i{A„(a;)tA„,(^),p}) 

is a dissipator and 

Mls = (a;)A„(a;)tA„Kw) 

a,Q:' 

is a Lamb shift Hamiltonian (commutes with Ms). 

(jj are Bohr frequencies, i.e. differences between eigenen- 
ergies of Ms- The summation is performed over all such 


( 2 ) 

(3) 

(4) 
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differences (positive and negative). Here [•, •] and {•, •} de¬ 
note a commutator and an anti-commutator correspond¬ 
ingly. Further, for an arbitrary operator T, 

T{u)= Y. Pie)TP{e'), (5) 

e' —e—uj 

where e and e' are eigenenergies of Hs (the summation 
is performed over all such pairs with the difference w), 
P{e) and P{e') are projectors onto the corresponding 
eigenspaces. Finally, 7aa'(w) and Saa'i>^) are correlation 
functions of the reservoir (defined by the operators Ba 
and the steady state of the reservoir); we will not need 
their exact expressions. If the reservoir is in a thermal 
state with the inverse temperature /3, the following rela¬ 
tion takes place: 

laa'i-Uj) = exp(-j3u)"/a'a(uj)- (6) 

Note that the stochastic limit technique [5] allows to de¬ 
rive not only a master equation for the reduced dynamics 
of the relevant system, but a stochastic unitary dynamics 
of the large system (the relevant system and the reservoir 
together). 

So, a master equation can be microscopically derived, 
but this requires a diagonalization of the system Hamilto¬ 
nian. This approach (and the master equation itself) are 
sometimes called global, because it considers the system 
as a whole. 

But the diagonalization of the system Hamiltonian can 
be a difficult task. Even if not, the master equation can be 
complicated and, hence, difficult for analysis. To simplify 
the task, suppose that the system Hamiltonian has the 
form 

Hs = Ho + i^V, (7) 

where the diagonalization of Hq is easy, and a dimension¬ 
less constant v is small. A particular case which we are 
interested in is a network of sites weakly coupled to each 
other and to a reservoir. In this case, Hq = Hj is a 
sum of free Hamiltonians of each site and V specifies their 
couplings; the diagonalization of Hq is often trivial. 

Of course, this situation is studied in the well-known 
framework of perturbation theory [5]. Suppose that we 
know eigenvalues } and the corresponding eigenstates 
{|ei°^)} of Hq. The eigenvalues and eigenstates of H can 
be represented by the series 

Sn = + ..., ( 8 ) 

|en) = \e^°'')+v\e^n^)+v^\e^n^) + (9) 

where the terms can be found explicitly. Suppose that 
the perturbation vV does not change the multiplicities of 
the zeroth-order eigenenergies (the situation when this is 
not the case will be discussed later). Then, as a zeroth- 
order approximation we can just put + 0(iy) and 


\e-n) = + 0{v). Accordingly, equations (H))-® are 

changed to 

^ = -i[Hq + vV + \^H^°Ip\+X^V^^\p) + 0{\^v), (10) 

{w)pA‘;°'> {w)'< 

w oc,a' 

-i{AW(u;)tAi°)(u;),p}), 

Pls '5'c<a'(w)A^°)(«;)1'A^°^(w), 

w Ct^Ct' 

where w = are differences between zeroth-order 

eigenenergies = E, and 

r(°)(w) = Y p^°\e)tp^°\e') 

E' — E—w 

for an arbitrary operator T. Here P^'^\E) are projectors 
onto the zeroth-order eigenspaces. 

This is exactly what is called the local approach to open 
quantum dynamics: the inter-site couplings does not affect 
the system-reservoir couplings. 

We see that the local approach can be justified on the 
base of the global one in case when the inter-site couplings 
are weak. Also note that both A and v are small parame¬ 
ters. By this reason, terms of orders v and A^ are included 
in (uni), but terms of order O^X^v) are not since this is a 
higher order of smallness. 

We assert that all discrepancies between the global and 
local approaches reported by various authors fall into 
at least one of two situations: either the inter-site cou¬ 
plings are not small or the analysed quantities are of or¬ 
der 0{X^v), which is beyond the local approach’s range of 
validity. 

Let us give an example of the second situation. In [5], 
the authors consider a simple network with two sites each 
of which is coupled to its own reservoir in a thermal state 
with different temperatures (see the next section), derive 
the expression for the steady-state heat flux based on the 
local approach, and show that, for some values of param¬ 
eters, there is a heat flux from the cold reservoir into the 
hot one, which violates the second law of thermodynamics. 
But the steady-state heat flux found in [2] is proportional 
to X^v'^ (in our notations). Since the local equation (ITUl) is 
written only up to terms of order A^, the heat flux should 
be set to zero in this order of smallness, and, hence, there 
is no violation of the second law of thermodynamics! 

If we want to derive a nontrivial expression for the 
steady-state heat flux of order X^v^ on the base of the 
local description, we should take the right-hand side of 
equation with the eigenenergies and their eigenstates 
expanded up to the second order in (I5|)-(in|). We are going 
to do this in the next section. 

Based on their analysis, the authors of [5] conclude that 
“the local approach is only valid for local observables such 
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as the population of each node [site], and is not valid for 
non-local observables describing energy fluxes”. Here we 
want to precise this statement: local approach is valid 
for observables of order A^, but not of order 0{X^i'). Of 
course, in the model considered in [5], where all non-local 
observables vanishes whenever v = 0, this is equivalent to 
the cited statement. 

Other words, phenomenological consideration of the lo¬ 
cal approach does not give its range of validity and leads 
to discrepancies with the microscopically justihed global 
approach and to violation of the second law of thermody¬ 
namics. In contrast, the viewpoint of perturbation theory 
naturally gives the range of validity of this approach and 
offers a way of obtaining corrections to it. 

An interesting question concerns the relation between 
the small parameters A and ly (e.g., whether A ^ i/ <C 1 
or ^ A <C 1, etc.). In principle, any relation is possible. 
But we should take into account that the corrections to 
the master equation ([5]) proportional to higher orders of A, 
in principle, can be derived. First steps in this direction 
have been done in [5] in the framework of stochastic limit 
mentioned above. If we drop the term of order O(A^), the 
expansion of the eigenenergies and eigenstates in ©-(ini) 
up to the nth order makes sense only if <C A. 

Physically, n can be interpreted as the ratio of a charac¬ 
teristic shift of energy levels to the characteristic difference 
between the non-perturbed levels. Also A^ can be inter¬ 
preted as a ratio of the time scale of relaxation of the 
reservoir to the time scale of evolution of the state of the 
system. 

Let us also mention that a perturbative treatment of 
inter-site couplings in the local description was also con¬ 
sidered in m- But this was done in a different way, which 
in some cases (e.g. degenerate eigenenergies or degenerate 
Bohr frequencies, i.e. if there are different pairs of energies 
with the same difference) produces indefinitely increasing 
(with time) terms. 

For recent considerations of the degenerate case, see 

[Ill- 

Two-Site model. — Consider a model of two sites A 
and B each of which is coupled to its own reservoir in a 
thermal state with different temperatures Th and Tc, ac¬ 
cordingly {Th > Tc). The corresponding inverse temper¬ 
atures are Ph and Pc- The sites are either two two-level 
systems or two harmonic oscillators. The Hamiltonian is 

H = Hq vV -\- Hh -f He XHah + XHbc, 

where Hq = EAa^a + Esb^b is a free Hamiltonian of 
two uncoupled sites, V = a^b + is their coupling, 
Hh and He are free Hamiltonians of the reservoirs, and 
Hah = (a -I- a'l') 0 Rh and Hbc = {b + b^) 0 Re spec¬ 
ify site-reservoir interactions. Here a (a^) and b (6^) are 
annihilation (creation) operators for the sites A and B 
correspondingly. They satisfy the relations 

aa^ + Sa^a =1, aa + 5aa = 0, 

bb'^+Sb^ = 1 , bb + Sbb = 0, 


where (5=1 for two-level systems and (5 = — 1 for harmonic 
oscillators. 

Our goal is to obtain a nontrivial expression for the 
steady-state heat flux from the hot reservoir based on 
the second-order perturbations in (El)-® and to ascer¬ 
tain that it is positive, i.e. does not violate the second law 
of thermodynamics. 

If the sites are harmonic oscillators, there is no need for 
perturbation theory, because the Bogoliubov transforma¬ 
tion [12] maps the system Hamiltonian Hs = Hq -|- lyV 
into a sum of two uncoupled oscillators (though each of 
them is coupled not to a single, but to both reservoirs). 
This considerably simplifies the task, and it is not so hard 
to derive and analyse the global master equation, which 
has been done in |2] . The case of two-level systems is more 
intricate for the global approach, and our “perturbed local 
approach” can be useful. So, we restrict our analysis to 
this case (5 = 1). 

We suppose that Ea > Eb- As shown in [2], this is a 
necessary condition for the violation of the second law of 
thermodynamics if we use the local approach. Also, the 
opposite case Ea < Eb can be trivially mapped into this 
one, because, actually, we will use that Pc > Ph only at 
the last moment; the case Ea = Eb is less difficult and 
was analysed in [3] in the global approach. Moreover, for 
simplicity let Ea — Eb = 1; the case of an arbitrary dif¬ 
ference Ea — Eb can be reduced to this one by a rescaling 
of the parameters of the system Hamiltonian. 

In the considered case, the eigenvalues and eigenvectors 
of the system Hamiltonian up to the second-order in v 
terms are: 


£oo — 0, |eoo) — |00), eu — Ea + Eb, \^ii) 
£01 =Eb- iy\ |eoi) = (^1 - |01) - |10), 

Eio = Ea + , |eio) = ^ |01) ■ 


111 ), 


Here |00), |01), |10), and |11) denote the both sites in the 
ground states, only site B in the excited state, only site A 
in the excited state, and both sites in the excited states, 
correspondingly. 

If we substitute these expressions to @-®, we will 
obtain the master equation 


— —i[Ho + vV, p] + X^Vh{p) + X^Vc{p), (11) 


where 


'Es{p) — ljs(EjspEjg r,{EjsEjs,p} 

j=A,B ^ 


+ exp{-PsUJj) (^F^^pEjs - , 


( 12 ) 


s = h,c. For shortness, we have adopted the notations 
7 js = 'Ys^ujj) (where 7 ^ is a correlation function of the 
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corresponding reservoir), Fjs = Fs{uij), where Fh = a+a\ 

F, = h + h\ 

LUA = Ell — Eoi = ElO ~ Eoo = Fa + ^ 

= Ell — Eio = Eol — Eqo = Eb — ■ 

It turns out that 

Fah = (1 ~ i'^)a — iy{a^a — aa^)b, 

Fbh = + v{a) a — aa))b^ 

Fac = — va{b% — bb^), 

Fbc = (1 ~ v^)b + va{b% — bb^). 

Also relation dG]) has been used in (EH). Finally, we 
have neglected the Lamb shift Hamiltonian: it leads to 
modifications of the parameters of the system Hamiltonian 
{Ea, Eb and the coefficient before V); this does not play 
a crucial role in our problem. Though, in other problems 
it can play a significant role [13] . 

The heat flux from the hot reservoir is given by the 
formula 

Jh = TiiVhip,){Ho + vV)], (13) 

where p* is a steady-state density operator. The substitu¬ 
tion of m into (IT^ yields the expression 

Jh = + Jt\a^a) + Ji^^b^b) + Ji^\x)), (14) 


where 

= lAh{l - 2i^^)[l + exp(-/3,iWA)] 

+ 7ScJ^^[l + exp{-PcOJB)] 

-|-7AcJ^^[l -l-exp(-/3ca;A)], 

= 7bc( 1 - 2i/^)[l -I- exp(-,5cWB)] 

+ + exp(-/3/jWB)] 

+ + exp{-/3hU!A)], 

= 7b,, 1/2 [1 -h exp(-/3,,WB)] 

+ 7A/t(l - J^^)[l + exp{-/3hUJA)] 

+ 7Aci^^[l + exp(-/3ca;A)] 

+ 7Bc(l - + exp(-/?cWB)], 

= 7bc[ 1 + exp(-/3cWB)] - jAhl^ + exp{-PhU}A)], 

= jAh exp(-/3,,u;A)(l - 2 i/ 2 ) + jbcb'^ exp{-l3cUJB) 

+ -fAcv'^ exp{-l3cUJA), 

= 7Bcexp(-,5cWB)(l - 2z/2) jbhb^ exp{-/3hU}B) 

+ ^Ahi^^ exp{-/3hUJA), 


where 

Jh'’ = lAhEA exp{-(3hioA) + IbhEbi^^ exp{-PhUJB) 

- 'yAhEBv'^ exp{-j3hUJB) 

= 7A/i[l + exp{-f}hiOA)][EA (l - - Ebv^)] 

= i^'^ilAhEA[l + exp{-Phi^A)] 

+7b/i£'b[1 + exp{-PhUJB)]) 

= jAhvEA[l + exp{-PhOOA)], 

X = a'fb + a&t (in our case V = X), and (T) stands for 
Tr[p*T] for an arbitrary operator T. 

The averages (a'I'a), (&^6), {X), and (Y) = i{a'^b — ab'^) 
satisfy the closed system of equations (follows from master 
equation EH)-dm) rewritten in the Heisenberg represen¬ 
tation) 

^{a^a) = -A2A/(““)(a^a) + —M<^‘^^\X) 

dr ' ^ ' 2 ' ' 

-v{Y)+\^c^^\ 

+ v{Y)+\^c^^\ 

— {X) = 
at 

- -b (y) -b 2X^vc^^\ 

— (Y) = 2v{a^ a) — 2v{b^b) 
at 

-(A)-yM(^^)(r), 


c(-^) = 7A,,exp(-/?,,WA) - 7Bcexp(-/3cWB)- 

Note that if we drop the terms of orders X^v and X^v^ in 
EH) and EH, we will obtain the corresponding equations 
from [5]. 

All left-hand sides in EH are zero since the averaging 
is performed over a steady state. So, we obtain a system 
of linear equations. The substitution of its solution into 
(fHl) yields the final result 


Jh = X^v'^ IAcEa 


FIbkEb 


1 - exp(-£’AA/3) 
1 -b exp{EAl3h) 
exp(£’BA/3) - 1 


exp{EBf3c) + 1 


(16) 


+ oiX^iy‘^), 


A/3 = Pc — Ph- Here we have dropped not only terms of 
order 0{X^v'^), but also terms of order O(A^), because the 
master equation itself is derived up to order A^. If we, 
however, hold terms of order O(A^), the expression will 
be much more cumbersome, but this will not affect the 
sign of Jh- Also, as we noticed at the end of the previous 
section, we assume that X (otherwise, the derivation 

of terms of order X^v^ makes no sense without corrections 
to the master equation of order A^, which are greater in 
this case). 

In the case of an arbitrary difference XE = Ea — Eb > 
0, the resulting formula EH is still valid with v substituted 
by v/XE. 

We see that, since A/3 > 0, the heat flux from the hot 
reservoir is always positive in accordance with the second 
law of thermodynamics. 
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Remark. — We assumed that the perturbation vV in 
© does not change the multiplicities of the eigenenergies. 
The opposite case is more tricky. If there is a removal 
of degeneracy, the distances between the new energy lev¬ 
els are small (formally, they are infinitesimal as z/ —>■ 0) 
and, hence, may be less than or comparable to the line 
broadening caused by the interaction with the reservoir. 
In this case, the quantum master equation approach or at 
least the secular approximation used in the microscopic 
derivation of the master equation is, generally speaking, 
not valid [T1I51I5] . So, equation (I2|), which we are based on, 
may be inadequate in this case. 

A possible solution is to go one step back from the 
master equation in the Gorini-Kossakowski-Sudarshan- 
Lindblad [T4l[T5] (GKSL) form ([2]) to the more general 
Redfield equation [nisiiis] and not to apply the secular 
approximation for pairs of Bohr frequencies with small 
differences. The resulting master equation will be not of 
the GKSL form, but it will reproduce the dynamics more 
adequately (the adequacy of the secular approximation is 
discussed in [T3l[T7]). It is interesting to mention that, in 
[TS] . the local approach was introduced exactly to avoid 
the secular approximation due to small energy differences, 
but, nevertheless, to stay in the equation of the GKSL 
form. 

However, two extreme cases can be discriminated. If 
the line broadening is much smaller than the distances 
between the levels from this new multiplet (A^ <C 1), 
we can still use the quantum matrix equation in form ([2) , 
i.e., with the full secular approximation. But if we restrict 
ourselves to the zeroth-order expansion of eigenenergies 
and eigenstates (0)-® in the master equation (local ap¬ 
proach), we should ignore the degeneracy of zeroth-order 
levels. 

On the contrary, if the line broadening is much larger 
than the width of this multiplet (difference between the 
maximal and minimal levels of it), i.e. u <§; <C 1, 
the removal of degeneracy should be ignored, and the full 
secular approximation can be applied. 

So, in both cases, we stay in an equation of the GKSL 
form ®. 

Conclusions. — The microscopically derived global 
master equation for quantum dynamics in open networks 
can be too difficult to construct and/or to analyse. The 
local master equation simplifies this task, but its range 
of validity can be too restrictive. Here we have proposed 
an intermediate approach, which uses stationary pertur¬ 
bation theory to obtain corrections to the local approach 
up to a desirable order of the inter-site coupling constant. 
This technique bridges the local and global approaches to 
construction of quantum master equations. 

* * * 
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